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Abstract. - We investigate the static shear viscosity on the sol side of the vulcanization 
transition within a minimal mesoscopic model for the Rouse-dynamics of a randomly crosslinked 
melt of phantom polymers. We derive an exact relation between the viscosity and the resistances 
measured in a corresponding random resistor network. This enables us to calculate the viscosity 
exactly for an ensemble of crosslinks without correlations. The viscosity diverges logarithmically 
as the critical point is approached. For a more realistic ensemble of crosslinks amenable to the 
scaling description of percolation, we prove the scaling relation k = <j> ~ P between the critical 
exponent k of the viscosity, the thermal exponent f3 associated with the gel fraction and the 
crossover exponent of a random resistor network. 



Introduction. - As the gelation or vulcanization transition is approached from the sol side, 
a melt or solution of polymers becomes increasingly more viscous, suggesting a divergence of 
the static shear viscosity rj at the critical point. The transition is commonly interpreted as a 
signature of a percolation transition ||J ||, ||: When the concentration c of crosslinks, which 
bind different polymers together, is increased to its critical value c cr it a macroscopic cluster 
of polymers is formed. Even for chemical gelation, i.e. permanent crosslinking, the measured 
values of the exponent k of the viscosity r\ ~ (c cr it — c)~ k lie in a broad range. Adam et al. [|] 
find values 0.6 < k < 0.9 for polycondensation without solvent and radical copolymerisation 
with solvent. Finite frequency measurements by Durand et al. [pi also yielded results within 
the above range. Silica gels and epoxy resins have been investigated by Martin et al . Jfj] , Q and 
by Colby et al. || |). The results are in the range 1.1 < k < 1.7. For a review see [icf . 

Using the percolation picture a number of heuristic proposals have been made as to how 
the exponent k is related to the critical exponents arising in the scaling description ]TT[ | of 
percolation. The most common proposal, k = 2v — (3, has firstly been given by de Gennes 
[ fl2"| within a Rouse approximation. Here v is the exponent governing the divergence of the 
correlation length and (3 is associated to the gel fraction. This result has been re-derived 
and supported with additional arguments i n p[ |l3| , 0, [| . On the other hand, the underlying 
assumptions have been questioned M. In ]14[ de Gennes hinted at an analogy between the 
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viscosity and the conductivity of a random mixture of superconductors and normal conductors, 
see |jj for details. The resulting conjecture is k — s, where s is the exponent ruling the 
divergence of the conductivity. In |D| Kertesz argued for s = v — (3/2 in analogy to the 
Alexander-Orbach conjecture |l6j. This is, at least, a reasonable approximation for s. The 
resulting value for k happens to be exactly the one sugg ested in || | for Zimm dynamics 
as opposed to Rouse dynamics. In fact, jD| |3) point out that the factor two by which their 
proposals for k for the two dynamics in question differ fits nicely to the two experimentally 
observed regimes. 

As can be seen most clearly from the arguments given in [^| , the basis for a representation 
of k in terms of the exponents of percolation theory is to relate the longest relaxation time 
of a cluster of macromolecules to the cluster's size. This relaxation time, of course, does not 
only depend on the assumed dynamical model but on the internal structure of the cluster 
as well. The internal structure of percolation clusters is by now believed to be characterized 
by an exponent independent of (3 and i>, the latter ruling the behaviour of the clusters on 
l arg e scales. This believe stems from the failure of the Alexander-Orbach conjecture p6j, see 
[y_7| and references therein. More quantitatively, the spectral dimension d s of the incipient 
spanning cluster introduced in []l6| is only approximately but not exactly equal to 4/3 for 
spatial dimensions d < 6. The spectral dimension encodes the fractal nature of the cluster's 
connectivity without giving reference to its spatial configuration. This can be seen, for example, 
from the fact that d s parametrizes the Lifshits tail D(E) ~ £;( d =>/ 2 ) _1 of the density of states 
D(E) of the discrete Laplacian within the cluster. Cates Jl8| relates the relaxation spectrum 
of a connected fractal cluster of n polymers to D(E) and finds that the static shear viscosity 
within Rouse approximation scales as nA 2 ' da ) _1 . 

Since the different heuristic arguments yield competing proposals even for virtually the 
same set of assumptions, the appropriate way is to discuss clear cut models. The purpose 
of this Letter is threefold: (i) Following the classical Kirkwood approach Jl^] we first 
calculate the static shear viscosity of a monodisperse sol of phantom monomer chains, which 
follow Rouse dynamics and are permanently crosslinked by Hookean springs chosen at random 
[ pp| . Identifying the crosslinks with electrical resistors soldered together by the polymers, 
the polymer network may be thought of as a random resistor network. This allows us to 
establish an exact correspondence between the viscosity and the resistance of a random resistor 
network. (ii) For the simplest distribution of crosslinks without any correlations, which 
amounts to a mean-field like model of percolation, this correspondence is employed to derive 
an exact expression for the averaged viscosity, implying a logarithmic divergence at the gelation 
transition. (iii) For general distributions of crosslinks, assumed to be amenable to the scaling 
description of percolation, the exact correspondence between the viscosity and the resistance 
is employed to establish the scaling relation 

k = 4>-p. (l) 

The exponent <f> was first introduced in the context of random resitor networks. It governs 
the growth of the resistance 7t(r) between two points on the incipient spanning cluster, which 
are a large spatial distance r apart: lZ(r) ~ r^l v (l7[ ^l] The exponent 4> is related 

to the spectral dimension, according to <j> = vdf[(2/d s ) — lj, where df = d — (3/v is the 
Hausdorff-Besicovitch dimension of the fractal. 

Dynamic Model. - We consider a system of N linear, identical, mono-disperse polymer 
chains, each consisting of L monomers. Monomer s on chain i is characterized by its time- 
dependent position vector Rt(i, s) (i = 1, . . . , N and s — 1, . . . , L) in d-dimensional space. 
We are interested in shear flows and impose an external velocity field v" xt (r,t). Here, Greek 
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indices indicate Cartesian co-ordinates, i.e. a — x, y, z, . . .. We will always use a flow field in 
the x-direction, increasing linearly with y, i.e. vjf xt (r,i) = n(t)y and u" xt (r,t) = for a =/= x. 
We employ the simplest, purely relaxational dynamics [ Jl8j | 

1 ri FT 

d t R?(i,s) = -- + v^ t (R«(i, a), t) + ff (i, «)• (2) 

In the course of time, the monomers relax to the state dH/dH = such that their velocity is 
equal to the externally imposed velocity field, i.e. dtRt(i,s) = u™ xt (Rt(i, s), i) . The relaxation 

process is disturbed by thermal noise £ with zero mean and covariance s) s') = 
2C _1 8a, p 5i,i' S SlS > S(t — t'). Here, the overbar indicates the average over the realizations of the 
Gaussian noise £. The relaxation constant is denoted by (, and we use energy units such that 
ksT =1. In (0) the Hamiltonian H := Hw + U is supposed to consists of two terms. The 
first one 

, N L-l 

H ™ ■= 2/2 E E ( R ^' s + !) - R ( 1 ' s )) 2 ( 3 ) 

guarantees the connectivity of each chain, the typical distance between monomers being given 
by the persistence length I > 0. The second one models M permanently formed crosslinks 
between randomly chosen pairs of monomers Q := {i e , s e ; i' e , s' e }%L±. To constrain the relative 
distance between two monomers, participating in a crosslink, we choose a harmonic potential 

, M 

U: =2^ E( R (^ S e)- R «> S e)) 2 > ( 4 ) 
e— 1 

whose strength is controlled by the parameter a > 0. For a — > hard crosslinks can 
be recovered pCfl . Since the Hamiltonian is quadratic in the monomer positions, it can 

be expressed in terms of a NL x NL random connectivity matrix T according to H =: 
W2« 2 ) £m'=i ££*=i R (*> s ) ■ T (h s; s') R(i', s>). 

Shear viscosity and resistor networks. - Our aim is the computation of the intrinsic shear 
stress o a p(t) as a function of the shear rate Following Jl9|], we express the shear stress 

in terms of the force per unit area, exerted by the polymers 



N 



N 



Here, pa stands for the polymer concentration and Fj*(i, s) := —dH/dR?(i, s) is the force on 
monomer (i, s). Since the dynamic equation (Q) is linear in the monomer positions Rt, it can 
be readily solved. Upon inserting this solution into (|^), we find the linear relation 

a xy (t)= j\rG(t-r) K (r), G{t) := ^ Tr((l - E ) expj- , (6) 

where Eq is the projector onto the null space of the connectivity matrix T. For a constant 
shear rate k, the intrinsic viscosity r\ is related to the stress tensor via r\ := <T X y/(KPo) such 
that 



r°° Cn 2 



1-E, 



(7) 
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is given by the trace of the Moore-Penrose inverse p3| of T, i.e. the inverse of T restricted to 
the subspace of non-zero eigenvalues. Note that we have made the dependence of 77 on the 
realization Q of the crosslinks explicit. 

Each crosslink realisation Q defines a random labelled graph, which can be decomposed 
into its maximal path- wise connected components or clusters, Q = UtLx-A/fc. The associated 
connectivity matrix is of block-diagonal form so that 

K N 

tW^^vWk). (8) 

fc=i 

Here, iV^ denotes the number of polymers in the cluster A4 . 

Let us identify a bond between two neighbouring monomers on the same polymer as a 
resistor of magnitude I 2 /a 2 and a crosslink between polymers as a resistor of magnitude 1. 
The resistance measured between any connected pair of vertices (£, s) and (i', s') will be denoted 
by K(i, s; i', s'). From (g) and |24[ Thm. F] it follows that 



^) = 4^V2 E K(i,s;i',s')- (9) 

k (i,s;i',s')eA/fc 

Together with (|J) this constitutes the announced connection between the viscosity and the 
resistances in a random resistor network. 

We expect an observable, like the viscosity, to be self-averaging in the macroscopic limit and 
compute the average (77) over all crosslink realisations in the macroscopic limit N — > 00, M — * 
00 with the concentration of crosslinks c := M/N being fixed. It will be advantageous to 
reorder the sum in (^) by summing first over all clusters consisting of given number n of 
polymers and subsequently over all "sizes" n. Thus we obtain for the averaged viscosity 

00 

(n) = J2 nr n ( V ) n , (10) 

where (if) n :— r^ViV -1 J2k=i $N k ,n vi-^k)) denotes the average of 77 over all clusters of size n 
and t„ := {N^ 1 X^feLi ^N k ,n) denotes the average number of clusters of size n per polymer. 

Completely random crosslinks. - The simplest distribution of crosslinks omits correlations 
between crosslinks and gives equal weight to all possible crosslink realisations. This case 
has been studied extensively in the theory of random graphs, as developed by |2q| . More 
precisely, we will calculate (//} := limjv- >00 (Ile= i( NL y 2 J2? c ,i> c =i 12s e , s ' e =i)v({ie, s e ] i' e , s' e }) in 
this section. In the macroscopic limit there are no clusters of macroscopic size for c < c cl -it := _^i 



and all monomers belong to tree clusters without loops [J25], Thms. 5d,e]. Moreover |25 
eq. (2.18)], one has t„ = 7i n ~ 2 {2c e~ 2c ) n /(2cn!). Since all L 2 (™- 1 ) n r! - 2 trees T n of size n 
can be proven to occur equally likely, the average with respect to them reduces to (r/}„ = 
L 2 ( 1 ~ n ^n 2 ~ n ^2 r T](T n ). Note further that the resistance between two vertices in a tree 
simplifies considerably, because there is a unique path connecting them, implying that all 
resistors are in series. In the simplest case of crosslinkcd Brownian particles, i.e. polymers 
consisting of just one monomer each (L = 1), the resistance 7Z(i; i') is equal to the number of 
crosslinks connecting i and i' . Hence, we refer to p6| , Thm. 1] for computing 

* — ' (n — v)\ 
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The double series arising from (|l0|), (Q) and (pd| ) can be summed up in closed form so that we 
obtain the averaged viscosity of crosslinked Brownian particles for all < c < | . This result is 
readily generalised to networks of crosslinked polymer chains (L > 1). The only difference is 
that now there are two different kinds of resistors with magnitudes I 2 /a 2 and 1, corresponding 
to intra-polymer bonds and crosslinks, respectively. Due to the fact that the monomer labels 
s are distributed independently from the chain labels i, one can reduce the averaged resistance 
to the Brownian-particle case (|lT|). This gives rise to the exact result 



Ca 2 f T I 2 n 



2 



lu| — !— ) 2c 

1 - 2c 



2 



— (L 2 -l) (12) 
12d y ' V ' 



for the averaged viscosity of a crosslinked polymer melt on the sol side < c < \ . It exhibits 
a logarithmic divergence at the critical concentration c cr it = \ corresponding to the critical 
exponent k = 0. The result (|l2|) is universal in the sense that the details of the model only 
affect the pre-factor of the critical divergence, which depends on the persistence length I, the 
"extension" of the crosslinks a and the length L of a polymer chain. For long polymers the 
viscosity is proportional to L 2 , as it should be within a Rouse- type model [^9|. This scaling is 
not altered when passing to the limit a J, of hard crosslinks. 

The disorder average of the viscosity has also been computed |27j] with the help of the 
replica trick for the simplest distribution without any correlations as well as for the crosslink 
distribution of Dcam and Edwards p8| which is believed to be more realistic. In both cases 
we find a logarithmic divergence of the viscosity for a network of Brownian particles, as the 
gelation transition is approached. 



Scaling description. - The above crosslink distribution ignores all spatial correlations and 
hence cannot account correctly for the critical behaviour of a three-dimensional crosslinked 
polymer melt. To improve on this shortcome we now assume that the mechanism of crosslink- 
ing generates clusters amenable to the scaling description of percolation. More specifically, 
we assume r„ = ?i~ T /((c cr i t — c)n <T ) with the scaling function / decaying faster than any 
polynomial for large arguments and being virtually constant for small arguments. The critical 
exponents a and r are related @, [ll) to the exponents (3 and v by a^ 1 = dv — f3 and t = 1+dva. 
We recall that the exponents /3, v and </> have been defined in the introduction, the latter ruling 
the behaviour of the resistance. 

It has turned out in the preceeding section that crosslinked Brownian particles and cross- 
linked polymer chains have identical critical behaviour. We adopt this simplification by 
assuming that it holds true in the present case, too. Using (||), we thus get 

\ fc=i i,ieA4 / 

The second assumption is that (77) , which only tests cluster of given finite size n, does not 
acquire any irregularity as the critical point is approached and that (f]) n \ c =c ■ ~ nb w ith 
some yet unknown exponent b. Note that, again, this holds true with b — 1/2 in the case of 
completely random crosslinks, see (|lT|). Furthermore, we remark that this assumption can be 
circumvented, at least to second order in e = 6 — d, by using the more sophisticated methods 
from p2]| . To determine b we use 

53 E ^ (v) n - ( E mj)) c ^ - cr-^-t (w) 

^ n=2 \j6Af(i) / 
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Here Af(i) denotes the cluster which contains particle i. In the first step we have used 
enumeration invariance. In the second step |l7], eq. (2.45)] has been employed. It follows 
that b — a<f> = (2/d s ) — 1 and (|l(]) finally leads to the scaling relation (|l]) for the averaged 
viscosity. 



Discussion. - The mean-field approximation |29[ ] for the percolation transition reproduces 
the critical behaviour |2j| of completely random crosslinks. The value k = for d > 6, which 
is found from (|l|) with the help of [jn], suggests a logarithmic divergence of the viscosity 
in the latter case. This is exactly what is found from the rigorous calculation. Moreover, 
the results of jp, || plugged into (Q) yield the e-expansion k = e/6 + lle 2 /1764 + C(e 3 ), 
e := 6 — d. For d = 3 we use the results of high precision simulations ]3(j in (|l|) and obtain 
^L— 3 ~ 0-71- In order to relate (Q) to the previously suggested scaling relations mentioned 
in the Introduction, we note that if (2/d s ) — 1 = (2/df) were true (which in general is not), 
one would recover de Gennes' original proposal k = 2v — (3. On the other hand, if the 
Alexander-Orbach conjecture, d s — 4/3, were inserted into 4>, the resulting scaling relation 
would be k = (dv — 3/5) /2. De Gennes' second porposal k = s together with Kertesz' result 
s = v—[3/2 yields yet a different scaling relation. This is even more surprising, because Kertesz' 
result is compatible with the Alexander-Orbach conjecture in d — 2. For comparison we list 
the values s\ , « 0.73, (2v — 0)\ ,_ 1.35 of the competing proposals for Rouse dynamics. 
The values are taken from [llj, table 2] . The similarity between the numerical values of k and 



s for d = 3 has to be interpreted as accidental, because for d — 2 one knows from duality 1 31 



that 4> = s, which yields with the high precision data from |32j the clearly distinct numerical 
values k\ d _ 2 ps 1.17 and s| 2 « 1.31. Thus the analogy of the viscosity to the conductivity 
of a mixture of normal and superconductors is incompatible with the interpretation advocated 
here. 

The above result k\ , ~ 0.71 agrees well with the experiments of [|J [5). On the other 



hand it is not compatible with the possibly oversimplifying albeit attractive proposal |2j, 13 
to interpret the wide variation of the exponent k of the viscosity as the signature of a splitting 
of the static universality class of gelation into different dynamic ones. Since Rouse and Zimm 
dynamics are considered |33|, to be at the extreme ends of the strength of the hydrodynamic 
interaction, the actual value of k should lie in the range bounded by the results for the 
respective dynamics. The divergence of the viscosity in the Zimm model has to be expected 
weaker than in the Rouse model, because the longest relaxation time in the Zimm model 
grows with a smaller power of the size of the cluster than in the Rouse model []l9|, |[ |) . Since 
our result for Rouse dynamics falls well into the range of small exponent values, we expect 
that the extension of our analysis to Zimm dynamics will yield a change for the worse, as 
far as the agreement with experiment is concerned. One may think of several other effects 
which may affect the predictions for k and have not been taken into account by our simple 
model. Excluded- volume interactions are known to be relevant even for thermostatic questions 
[ p3[ [llj. We refer to |l^, || and references therein where attempts are made to account for 
these effects on a heuristic level. Additionally, one expects entanglement effects to play a vital 
role in stress relaxation p9] . However, the static viscosity measures stress relaxation only on 
the longest time scales in the sol phase. It has been argued || |34| that in the regime close to 
the transition entanglement effects are less important because of two reasons: first, there are 
almost no permanent entanglements in the sense of interlocking loops. Second, the time scale 
of a temporary entanglement of two clusters is determined by the smaller cluster, whereas the 
dynamics on the longest time scales is determined by the larger cluster. 

We hope that the analysis of the simple model presented in this Letter will serve as a 
first step towards a more sound understanding of dynamical effects in gelling polymeric sols. 
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Clearly, further investigations are necessary to embed the variety of well-developed scaling 
pictures within an appealing theoretical framework. 
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